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We reconsider both the global and local stability of solutions of continuously evolving 
00 ' dynamical systems from a geometric perspective. We clarify that an unambiguous definition 

o ■ 

' of stability generally requires the choice of additional geometric structure that is not intrinsic 



to the dynamical system itself. While global Lyapunov stability is based on the choice of 
seminorms on the vector bundle of perturbations, we propose a definition of local stability 
based on the choice of a linear connection. We show how this definition reproduces known 

oo 

' stability criteria for second order dynamical systems. In contrast to the general case, the 

special geometry of Lagrangian systems provides completely intrinsic notions of global and 
O |. local stability. We demonstrate that these do not suffer from the limitations occurring in 

■ the analysis of the Maupertuis- Jacobi geodesies associated to natural Lagrangian systems. 
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\0 ; I. INTRODUCTION 



Continuously evolving dynamical systems universally appear as models in every branch of nat- 



OO ! ural science. Their usefulness is determined by their predictivity which in turn hinges on the 

o 

stability of their solutions. Given some uncertainty in the measured initial conditions, a physically 
^ , meaningful system should offer control on the possible evolution of deviations from a given refer- 

' ence trajectory. Global control of late-time deviations is as important as a local understanding of 

stability. An extensive mathematical theory exists for the global stability of solutions of dynami- 
cal systems; the fundamental quantities in this theory are the Lyapunov exponents that measure 
exponential deviations . These exponents are often very hard to determine analytically, and 



Sometimes 



various approaches for their calculation are used in the literature 
the stability of the same dynamical system is considered from different points of view, and the 
corresponding Lyapunov exponents are compared; this leads to some debate on the question of 



which point of view is preferable 
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141 ] . The local stability of solutions of dynamical 
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systems on the other hand seems less well understood. To our knowledge, the existing mathemat- 
ical theory is applicable only to second order differential systems. It was developed by Kosambi, 



Cartan and Chern 



17l | , and covers as a special case the local stability analysis of geodesies 



on a Riemannian manifold via the Jacobi equation. 

In this article we describe dynamical systems and their perturbations in a unified geometric 
language which is reviewed in section [Til This language has the advantage to specify very precisely 
the intrinsic geometric ingredients that define the dynamical system. It is worth emphasizing that 
the translation of dynamical systems into geometry is unambiguous. However, this is not true 
for the definition of stability. We will clarify in section IIIII that stability in general is not an 
intrinsic notion but requires geometric structure not included in the definition of the dynamical 
system. This fact is not sufficiently appreciated in the literature. We will see that global Lyapunov 
stability entails the choice of seminorms on the vector bundle of perturbations. We will then 
propose a general definition of local stability applicable to any dynamical system. This notion of 
local stability requires the choice of a linear connection, and reduces to standard stability criteria for 
second order systems. In contrast to the general case, Lagrangian systems offer intrinsic geometric 
structures to define global and local stability, as we will discuss in section IIVI We will then consider 
natural Lagrangian systems in section |V] for which independent notions of stability exist based on 
the translation of solutions into geodesies via the Maupertuis-Jacobi principle. We point out that 
this detour suffers from a number of limitations which are nicely resolved in the direct intrinsic 
approach valid for all Lagrangian systems. We conclude in section IVTl 



II. PERTURBATIONS OF DYNAMICAL SYSTEMS 



In this section we review the basic geometric structure common to all dynamical systems and 
their perturbations. The geometric language used here will form the basis for our discussion of the 
stability of solutions in the following sections. 

Fundamentally, deterministic dynamical systems can be introduced as formal rules that describe 
the evolution of points in some set S with respect to an external, discrete or continuous time 
parameter running in another set T. More precisely, a dynamical system is a map 

X S ^ S , {t,x) ^ <)){t,x) (1) 

which satisfies </'(t, •) o '/'(•s,-) = 4>{t + s, •) for all times t,s G T. This bare definition must be 
enriched with additional structure if we want to model realistic dynamical systems. 
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In the following we will therefore concentrate on the important class of continuously evolving 
dynamical systems for which T = M. Geometrically, we focus on point sets S that are differentiable 
manifolds, and assume that a vector field X over S is given: 

X :S ^TS,p^XpeTpS. (2) 

The integral curves 7p of the vector field X through a point p £ S are defined by two conditions: 
that their tangent vector 7p(t) = X(7p(t)) agrees with the vector field at each point, and that 
7p(0) = p. For simplicity we assume that these integral curves are defined for all times. Then a 
dynamical system in the sense of ([1]) is defined by the fiow of the vector field X, i.e., by setting 
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FIG. 1: Trajectories of dynamical systems and their perturbations. 

Any integral curve 7p of the flow vector field X provides a solution trajectory of the dynamical 
system with initial condition p. Perturbations of 7p are studied by varying this initial condition. 
Any curve p{(j) with j>(0) = p provides a continuous family of initial conditions which evolve under 
the dynamical system as (j){t,p{a)) = 7p(o-)(^)- All curves so obtained are continuous variations 
of function of the two parameters t and a we thus obtain a two-dimensional surface 

(j){t,p{a)) in S, see figure [TJ In consequence the two tangent vectors X and i{t,a) = 9(j|<^(t,p{o-)) 
commute, [X, ^] = 0. Along the reference solution 7^ which lies in this surface at o" = we 
can now study the evolution of the tangent vector ^(t) = d(j\^^(^t)- Precisely this vector governs 
the linear stability of the solution 7^, i.e., its stability with respect to infinitesimal perturbations. 
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Globally the perturbation vector field ^(t) along 7^ is determined as the push-forward of an initial 
perturbation vector ^0 £ TpS via 

m = m,-)Uo. (3) 

To summarize, the dynamical systems (S, X) that we will study in this paper have a very simple 
geometric structure: 

• they are defined in terms of a flow vector field X over some differentiable manifold S; 

• their solutions are integral curves of X; and 

• the linear perturbations of a solution 7p are vector fields ^ along 7p for which [X, ^] = 0. 

Choosing local coordinates on (S, the dynamical system thus becomes a system of first-order 
differential equations over S for the integral curves of the vector field X, 

x(ty = X\x{t)) . (4) 

Expanding the commutator condition [X, ^] = for the linear perturbation gives 

dX 



e{t) 



e{t). (5) 

x{t) 



These local expressions seem to suggest that the above framework for dynamical systems would 
only allow for first order differential equations on a configuration space. However, this crucially 
depends on the choice of manifold S, its relation to the configuration space, and the type of flow 
vector field X. To see this we now discuss in the same geometric framework the case of second 
order dynamical systems which will play a major role in this paper (higher than second order 
construction require higher order tangent or jet bundles, see jisl). 

Hence consider a system whose configuration space is modeled by a smooth manifold M. Initial 
conditions for a second order system consist of an initial position q £ M and an initial velocity 
V G TqM. Combining these into pairs {q,v), one realizes that a natural geometric arena for second 
order systems is the tangent bundle TM. It is clear that not all vector fields over TM can produce 
second order dynamics on M. Indeed, a second order system is defined by a vector field X on 
S = TM with the additional condition 



TT^X = u 



(6) 
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for all (g, u) S TM, where tt* is the push-forward of the projection map vr : TM — > M , (g, n) h-> q. 
This condition defines X as a so-called semispray; in local coordinates {x'^,u'^) on TM one finds 
X = u""-^ + ^2 afs:- This ensures that the integral curves of X, which satisfy 

x'^(t) = n'^(t) , u''{t) = X^{x{t),u{t)) x^{t)=X^{x{t),x{t)), (7) 



represent a second order differential evolution on the configuration space M. In geometrical terms, 
the integral curves 7 of X become natural lifts 7(t) = (7(i),7(i)) to TM of curves 7 in M. 

Once second order systems on M have been rewritten as (special cases of) first order systems 
on S" = TM, perturbations of solutions can be analyzed as before. The perturbation 

along a given solution 7 will now be a vector field over TM with [X, ^] =0. It is instructive to 
expand this condition in local coordinates on TM, which yields 

ca — ca CO. — S> , cb , i'a _ cb , rQ\ 

^1-^2, ^2-^?l + ^?2 => ^1 - + • 

As expected, the perturbation equation can be read as a second order differential equation. Hence 
we see that the simple geometric framework for dynamical systems and perturbations of solutions 
can be easily applied to higher order dynamics. 



III. GEOMETRIC NOTIONS OF STABILITY 



In the previous section we have defined dynamical systems (S, X) via vector fields X and their 
associated flow on some differentiable manifold 5. We have also presented the geometric language 
to define linear perturbations of specific solutions as connecting vector fields that commute with X. 
In this section we discuss the issue of the stability of solutions under such linear perturbations, both 
globally and locally. We will find that stability generally is not intrinsic to a dynamical system, 
but requires additional geometric structure beyond S and X to be well-defined. 



A. Global stability 



Consider a solution 'jp^ of a dynamical system {S,X). The linear perturbations of such a 
solution are described by tangent vector fields ^ along 7 and satisfy [X, ^] = 0. Given an initial 
perturbation ^0 £ ^po*^' dynamical system completely determines the evolution of along the 
solution 7 via the push-forward as in equation ([3]). Since, intuitively, ^ points to a neighboring 
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solution, we wish to call 7 unstable, if any choice of .^o results in exponential deviations of ^ for 
late times. 

In order to implement this idea, we proceed in two steps. First we need to quantify the size of 
a perturbation ^; so the tangent bundle TS must be equipped with a family of (semi-)norms 

||-||p:rp5^]R+,e^||e||p (10) 

for each point p of 5. This is additional structure not provided in the definition of the dynamical 
system {S,X). Second, once this geometric structure is available, we are in the position to define 



the so-called Lyapunov exponent [1 



A(po,Co) = limsup jln||(0(t, •))*ColU(t,po) = limsupyln||^(t)||^ (f) (11) 

for given initial point po € S and initial perturbation vector ^0 £ '^po^- Exponential growth of a 
perturbation will result in positive Lyapunov exponent A; therefore, a given trajectory identified 
by pq is defined as stable if and only if X{po,£,o) < is non-positive for all initial perturbations .^o- 
This is a definition of global stability because the Lyapunov exponents are evaluated for late times. 

We note that the Lyapunov exponents and the notion of global Lyapunov stability can be 
straightforwardly applied to any dynamical system formulated as a vector flow. For instance, in 
the case of a second order system over a configuration space M, which we have seen is a vector 
flow over S = TM, the additional geometric structure needed to define global stability is a family 
of seminorms on the bundle TS = TTM. We will see this construction at work below. 

Even though a family of seminorms is sufficient to define (jlip , many general results in the theory 
of Lyapunov exponents have been achieved by equipping S with a Riemannian metric and consid- 
ering the associated family of norms on TS. In this case, it can be shown that for a given point p 
the Lyapunov exponent \{p,C) can assume at most n distinct values, < . . . < Moreover, 
there exists a so called filtration Ei <Z . . . d Es = TpS such that X{p,C) = ^ov ^ Ei\ Ei-\. 
So initial perturbations from this subset of TpS result at late times in exponential deviations con- 
trolled by the corresponding Lyapunov exponent. Stability theory based on Lyapunov exponents 
is a vast subject; we will not detail its major achievements here but refer the interested reader 
e.g. to [3|, yj. We only remark the following notable result: on a compact Riemannian manifold 5, 
the supremum limit in (jlip can be replaced by the ordinary limit which, moreover, is finite and 
independent of the particular Riemannian metric defined on S. 

One could suspect that, even though the Lyapunov exponent quantitatively depends on the 
specific choice of seminorms on TS, the qualitative notion of stability would not. The actual 
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independence of Lyapunov exponents from a Riemannian metric on compact manifolds seems to 
favor this hypothesis. However, the following simple example immediately reveals that this is not 
the case. Consider the inverted harmonic oscillator over M = M with Lagrangian 

L = i(u2 + ^V) (12) 

for fi > 0. It is easy to see that the associated flow vector field on TM is given hy X = {u, fi'^x). A 
generic perturbation of the solution x = Ae^^ + Be~^^ is given by ^i(t) = (ae^* + be~^^)dx, where a 
and b are arbitrary constants parametrizing the possible choices for the initial perturbation vector 
^ G TpTM. The point p at t = has coordinates (x,u) = (A + B,^(A — B)). The simple choice 
of seminorms induced on TM by the Euclidean metric g = dx®dx gives the Lyapunov exponents 

A(A, B, a, b) = lim - ln(ae'^* + fte"'^*) . (13) 

t-^oo t 

The Lyapunov spectrum is then given by {—fi,fj,}, where —fi is associated to the subspace a = 
of TpTM and ^ to its complement. Therefore the solution is classified as unstable. However, why 
should we not consider the metric g = (x^ + \)~^dx ® dx? Then the Lyapunov exponent 

1 rtpM^ _|_ hp-i^^ 

\{A, B, a, b) = lim - In , ^ ^ (14) 

vanishes (excluding the fixed point case A = B = 0). The solution is now classified as stable in 
contradiction to the previous result. 

We emphasize that the key point behind the Lyapunov notion of global stability of the trajec- 
tories of a dynamical system {S, X) is the necessity for additional geometric structure: seminorms 
on TS that can measure the size of the perturbation. 



B. Local stability 

We now turn to a very intuitive proposal for the local stability of solutions of any dynamical 
system {S,X). We wish to classify a solution as locally stable if the acceleration between ^p^ 
and a neighboring solution is negative, i.e., if the curves are forced together. To implement this 
technically, we need to be able to calculate the second derivative of the perturbation vector field ^ 
that points to the neighboring solutions. The only derivative that is available without introducing 
further structure is the Lie derivative, but already the first derivative CxS, = [X,S,] = vanishes. 
So we need to require another derivative. 

The definition of local stability again proceeds in two steps. First, we introduce a linear connec- 
tion VxY that sends a pair of vector fields X, Y over S to another vector field over S. Using the 
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connection, one finds that the second covariant derivative of ^ along the flow X is a linear operator 
on TS which can be written as 



yx{T{x,-) + y.x) + {T{x,-) + v.xf ■i = R{x)-i, (15) 



where T{X, Y) = — VyX — [X, Y] is the torsion tensor of V, and both T(X, •) and V.X are 

linear operators on TS. Second, we define a solution of the dynamical system as locally stable 
if all eigenvalues of the linear operator R evaluated along are non-positive. 

To our knowledge this definition of local stability, which can be applied to dynamical systems of 
any order, has not appeared in the literature. However, we will now show that our proposal has a 
very nice relation to KCC stability theory as developed by Kosambi, Cartan and Chern 15l.ll6l. Il7[ 



for second order dynamical systems over a configuration space manifold M. We have already seen 
that these systems are defined as vector flows over the tangent bundle S = TM in terms of a 
semispray X. Without loss of generality we may use local coordinates = {x"',u°') on S to write 

X = u''^-2G^{x,u)^. (16) 

To discuss the geometry of the dynamical system so defined, KCC introduce a non-linear con- 
nection on M with coefficients A^^^ = This can be understood in terms of a dynamical 
covariant derivative V^: for two vector fields v,w over M, 



v''da{w'') + N\{x,v)w^]-^^. (17) 



(In the special case of a linear connection one simply has N"'i,{x, v) = T"' i)c{x)v'^ .) Using one can 
rewrite the second order equation in Q, which governs the evolution of a perturbation vector ^1 
over M, in the form 



V,^Vf 6 = Pix) • ei (18) 



for a linear operator P on TM the components of which may be found e.g. in jl9(]: 

A solution of the dynamical system is called KCC-stable if all eigenvalues of P evaluated along 
7po are non-positive. 

Now to define local stability of the system {S = TM, X) with X provided by (|16p according 
to our proposal, we need to specify a linear connection V over S = TM instead of a non-linear 
connection over M. For this purpose we note that for every non-linear connection A^ on M 
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exists an associated compatible linear connection V on TM, which is said to be of Berwald type. 
Technical details and the global form of this linear connection are given in [2^ . The components 
of V are defined by V aEb = ^^a-^c with respect to a basis (Ea) of TTM. Choosing the basis 
(Ea) = (gfs- — N^a^jb^ adapted to A^, the only non-vanishing connection components are 

r = r^- = f20) 

In the same basis, we now calculate the linear operator R{X) on TS = TTM defined in (jl5p . the 
eigenvalues of which are relevant for our proposal for local stability. The result has the form 



^2 b b + ^ b 



(21) 



in terms of the deviation tensor P that appears in the stability equation (jlSp of KCC theory, of 
e"" = 2G'^ — N'^bU^, of {deYh = and of two more tensors Ri and R2- We only need to know that 
9e = implies Ri = 0; in this case the eigenvalues of R are precisely those of P, so our definition 
of local stability reduces to the KCC definition. Note that (9e = is solved by 

G»(x,«) = G^2)(x,n) + ie'^(x) (22) 

where G^2)(2;,u) is homogeneous in u of degree two. We will see that semisprays with such G°' 
arise from natural Lagrangian systems. So our criterion for local stability agrees with the KCC 
criterion for all natural Lagrangian systems. 

The fact that local stability of a second order system can be equivalently understood in terms of 
a nonlinear connection on M or in terms of a linear connection on TM as we propose was already 
shown for e = in [2^, but in a different way. This is in fact the only case where the KCC stability 
treatment seems to be natural, because the second order system then takes the autoparallel form 

i: = 0. But also the Berwald type connection satisfies = 0. For all other cases we do not 

consider a non-linear connection a particularly natural structure. Another point in favor of our 
linear connection proposal is the fact that it is naturally defined for any dynamical system (5, X) , 
not only for those related to second order dynamics. 

We emphasize that the key point behind the notion of local stability of the trajectories of 
a dynamical system (S*, X) here proposed is the necessity for additional geometric structure: a 
connection V on TS that allows to measure the acceleration between nearby curves. 



10 



C. Geodesies 

Both the definitions of global and local stability of the solutions of generic dynamical sys- 
tems (S, X) require the choice of additional geometric structure: in the first case seminorms have 
to be chosen on TS to measure the size of the perturbation; in the second case a connection on TS 
has to be chosen to measure local accelerations. 

However, some dynamical systems may provide natural choices. In the case of geodesies on a 
Riemannian manifold (M, g), which can be seen as a second order dynamical system (S = TM, X), 
natural choices exist for the seminorms as well as for a connection on TS = TTM. 

The seminorms for the definition of global Lyapunov stability can be constructed from any one 



of the well-known lifts of the metric g to the tangent bundle which are discussed in detail in 
As an example we may choose seminorms 

iieii^ = (5^)p(e,e) (23) 

on TTM by employing the vertical lift of g to TM which is defined by g^ {X,Y) = (/(vr* X, vr^y) 
for vector fields X, Y on TM. This construction relates instability to the exponential divergence 
of trajectories purely in the configuration space. In the notation of ([8|) the degenerate bilinear 
form g^ projects out the components Indeed, for a geodesic 7 : M — > M with initial condition 
(g, v) = p ^ TM and initial perturbation (^1, ^2) = ^ G TpTM, the Lyapunov exponents are simply 
calculated as 

A(p,0 = lim lln5^(t)(6(t),6(t)) • (24) 

Affinely parametrized geodesies with tangent vector field x over M obey the equation V|''^i; = 
where V^*^ is the Levi-Civita linear connection over M determined by the metric g. As discussed 
in the previous section, there is a linear connection V over TM naturally associated to V^'-'' (which 
now plays the role of N). The connection V of Berwald type is nothing else but the horizontal 
lift V = of V^^ to TM; see 2l|] for details. According to our definition, local stability of a 
geodesic requires that the eigenvalues of R{X) as defined in (jlSp should be non-positive. From 
the calculation made above in the comparison to KCC theory we hence deduce that local stability 
becomes equivalent to the eigenvalues of P{x)"'h = —R^pbqX^x'^ being non-positive. Considering 
the geodesic deviation equation 

(Vi^Vi^)'^ = -R%tgxPx'^ei (25) 
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for a simple perturbation [x, ^i] = on M then tells us that our notion for the local stability of 
geodesies naturally reduces to the standard definition of geodesic stability. 

So geodesies admit natural definitions of global as well as local stability. In the following section 
we will show that also Lagrangian systems offer intrinsic notions of global and local stability. 



IV. STABILITY ANALYSIS OF LAGRANGIAN SYSTEMS 

In the previous sections we have introduced all necessary ingredients for a discussion of the 
stability of solutions of Lagrangian systems on some configuration space manifold M. In order 
to apply the geometric machinery we will first formulate these systems as vector flows over TM. 
Then we will demonstrate the existence of natural notions of stability that do not require further 
information apart from that given by the Lagrangian system itself. 



A. Geometric formulation 



A Lagrangian system is a pair (M, L) of a configuration space manifold M and a Lagrangian 
function L : TM R. Locally this is a function L{x, u) of coordinates on M and coordinates 
on TxM that give the components of a vector u in the coordinate- induced basis The trajectories 
of the Lagrangian system are the curves 7 : M — > M with tangent vectors 7 for which the variation 
of the action integral vanishes, 

= 55[7]=5 r L{j{t),j{t))dt. (26) 
Jti 

In the variation the curve's endpoints 7(^1) = pi and 7(^2) = P2 are held fixed. Once coordinates 
have been chosen, solutions x"'{t) must satisfy the Euler-Lagrange equations 

d dL dL . , 

dtdx"- dx°- ' ^ ' 

As any second order dynamical system, a Lagrangian system can be described as a vector 
flow X on S = TM. Since this formulation is not so well-known, we present the relevant steps 
following 



22 1 . First, one defines the Liouville vector field A = u'^^f-^ on TM. This produces the 



du 

energy function E : TM M hy E = A(L) — L. Now recall that a vector field Y on TM is said 
to be vertical if it is tangent to the fiber, n^Y = 0. Moreover, a vector w = w"'-^ on M, has 



an associated vertical lift = w'^-^ to TM. Combining the vertical lift and vr^, defines a (1, 1 

du- 



tensor field S on TM by S{Y) = (vr^y)^; in coordinates, S = ® dx"'. With these ingredients. 
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we can construct the two-form to = d{dLo S) on TM. As a final step, the vector field X generating 
the Lagrangian flow is the unique vector field X solving 

ixuj = -dE . (28) 

This construction is easily understood in a local chart (x, u) of TM. Simple algebra shows that 

dE = - — — rn" - ^— dx" + - — —ru^du'' , 
\ox"'ou'' ox"' J ou^-ou" 

oj = - — dx''^dx^ + - — —rdu^Adx^. (29) 
Writing X = Xf-£^ + X^-£^, condition ([28]) is equivalent to 



This shows that the vector field X on TM is a semispray, as needed for a second order dynamical 
system. The integral curves of X indeed solve the Euler-Lagrange equations. 

B. Intrinsic notions of stability 

Within the geometric picture of Lagrangian systems (M, L) as dynamical systems {S = TM, X) 
we are now in the position to ask for the stability of solutions. We will discuss global and local 
stability in turn. 

The definition of global stability via Lyapunov exponents requires the specification of a family 
of seminorms on the tangent bundle TS. A Lagrangian system provides us with an intrinsic object 
that can be used for this purpose, see e.g. [23]: the generalized Lagrange metric 



Although the coefficients of g depend on the fiber coordinates u for generic Lagrangians L, the 
generalized Lagrange metric transforms under coordinate transformations of 5 = TM as a (0, 2)- 
tensor on M. Regularity of the Lagrangian corresponds to rank g = dimM. The situation for the 
Lagrangian system is now very similar to that of geodesies on a Riemannian manifold. Instead of 
lifting the Riemannian metric to TM, we simply need to extend the generalized Lagrange metric g 
to TM. As for lifts, there are several possibilities. As an example we consider the case, where g 
as displayed above is simply read as a tensor over TM. For a reference solution 7 : M — > M with 
initial condition {q,v) = p £ TM and initial perturbation {^1,^,2) = C ^ TpTM, the Lyapunov 
exponents are then simply calculated as 

A(p,0 =/_im lln5^(,)(ei(t),6(t)). (32) 
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This is of the same form as the formula for geodesies where we used the vertical lift of the Rie- 
mannian metric. So this choice of seminorms measures the exponential deviation of perturbations 
purely within the configuration space M. 

To define the local stability of solutions of (M, L) we need to specify a linear connection acting 
on vector fields on 5 = TM . Again the Lagrangian system provides us with an intrinsic choice. 
We may simply use the Berwald type linear connection V on TM which is associated to any given 
semispray. The relevant quantities G"^ and A^"^;, = needed to determine this connection follow 
by combining (jl6p . (j30p and the definition of the generalized Lagrange metric: 



2" dvPdx" 2'' dx^' 



In the basis {Ea) = —N^a-ffr, gfs:) of TS the connection V is then provided by the connection 



coefficients (pOl) . 

Both the global and local stability of solutions of a Lagrangian system (M, L) can be defined 
intrinsically, without the introduction of further geometric structure. This makes essential use of 
the fact that the Lagrangian system can be rewritten in the dynamical system language as a vector 
flow X over S = TM. 



V. NATURAL LAGRANGIAN SYSTEMS 



The solutions of natural Lagrangian systems can be translated using the Maupertuis-Jacobi 
principle to geodesies on an associated Riemannian manifold. We have seen that both global and 
local stability can be easily defined for geodesies, and this fact is often used in the literature to 
analyze the stability of natural Lagrangian systems. The main purpose of this section is to clarify 
the differences between the intrinsic notions of stability discussed in the previous section and the 
notions of stability obtained from the translation principle. It will turn out that the intrinsic 
notions of stability have considerable advantages. 



A. Intrinsic stability 

Before we translate natural Lagrangian systems into geodesies, we specialize the intrinsic notions 
of stability developed for general Lagrangian systems to this case. A natural Lagrangian system is 
defined by a Lagrangian function L : TM — > M which is quadratic in the velocities: 

L(x,n) = ^fc,fe(x)nV-y(x). (34) 
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This corresponds to a dynamical system {S = TM, X) for which the flow vector field X over TM 
is determined according to the general expression (I30p . and in terms of the Levi-Civita connection 
of the kinetic energy metric k, as 

X? = , X^ = -IC" = -T^.fbcu'u' - V^kfV . (35) 

The intrinsic notion of global stability is based on the generalized Lagrange metric g of (j3ip 
which here agrees with the kinetic energy metric, g = k. Assuming a regular Lagrangian and 
positive energy amounts to the requirement of positive definiteness of k. The second order evolution 
equation of perturbations .^i can be rewritten as 

(v(^-),vw, Ci)" = - (%) V.i^i'^ + V(fc)«v('=),y)e? . (36) 

The calculation of Lyapunov exponents with respect to some extension of k to TM is unproblematic 
and follows ([32]) . 

The intrinsic notion of local stability is linked to the Berwald type connection on S = TM. 
Note that the natural Lagrangian system has a flow vector field X2 = —20°" of the type (j22p . so our 
general calculation above shows that local stability becomes equivalent to the KCC definition. In 
practice, therefore, local stability means non-positive eigenvalues of the endomorphism P defined 
in ([19]). Specializing to the case here we find 

P\ = -R^kfpb.uPu'^ - Vikf^^^'hV , (37) 

where u has to be evaluated for the tangent vectors x along a reference solution. This is the same 
operator as the one we found governing the evolution of perturbations. 

Reconsider the example of the inverted harmonic oscillator of equation (jl2p . The generalized 
Lagrange metric now gives a preferred intrinsic choice of seminorms in terms oi g = dx dx, thus 
leading to the expected global Lyapunov instability of the system. For local stability we use the 
fact that our stability criterion and the KCC criterion agree, since we consider a natural Lagrangian 
system. From ()37p we find P = fi^ > 0. This shows the local instability of all solutions under this 
intrinsic choice of connection. 

B. Riemannian translation 

Geodesies on a Riemannian manifold with fixed affine parametrization can be viewed as a 
Lagrangian system. The converse is also true for natural Lagrangian systems. We now state a 
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version of the Jacobi-Maupertuis 



e.g. in scalar field cosmology 
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principle, including indefinite kinetic energy metrics which occur 



25( 1 ; then we will use this translation principle to investigate how 
the geodesic stability and the intrinsic stability of natural Lagrangian systems compare. 

Let (M, L) be a natural Lagrangian system (j34p with non-degenerate quadratic form k. Then 
the solution trajectories of fixed energy 

E = habix)x''x'> + Vix) (38) 

of the Lagrangian system coincide with the geodesies of the metric manifold (M, g) with 

gab = C\E-V{x)\kab (39) 

for constant C > (the conventional choice is C = 2). The time parameter t of the Lagrangian 
system is related to the afRne parameter r of the geodesies, which we define by \gab^§p^§f\ = 1? as 



dt 



dr 



^ (40) 



V2C\E-V{x)\ 

To see why the Jacobi-Maupertuis principle holds, we write the Euler-Lagrange equations of 
motion of the Lagrangian system as 

x"" + r(^kfbcx''x^ + k''^dbV = 0. (41) 

Now consider geodesies on a metric manifold (M, g) with the conformally rescaled metric 
9ab = f^^ab) parametrized in terms of an arbitrary non-affine parameter also called t. These 
geodesies satisfy the equation 

x" + r^g)\,x''x' - h'^j^ In \gpgxPx'^\ = , (42) 

now using the Christoffel symbols T^g) of the rescaled metric g, and may be rewritten in terms of 
the original kinetic energy metric k and the conformal factor a as 

x-^ + r^„f,,x-x' - {kpgxPx'^) k-'db In \a\ + In = . (43) 



This equation is equivalent to the Euler-Lagrange equations (j4ip provided that 

k^^dbV + (k^ai^i") k^^db In lal - ±'^4 In ^ = . (44) 

dt \kpqxVxi\^/^ ^ ' 

The contraction of this condition with i'^/cca is easily integrated with integration constant 
yielding the Hamiltonian constraint ()38p . We solve this for k^qXpifl and substitute back into the 
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above condition for the equivalence of the geodesic and Euler-Lagrange equations. Rearranging 
the terms we obtain 

This is satisfied for all possible trajectories only if a'^ = C \E — V\ for positive constant C, which 
demonstrates the first part of the claim. Finally, \gab^^%f\ = 1 defines an affine parameter r 
along the geodesies of {M,g). We may now use the Hamiltonian constraint and our result for cr^ 
to show that (dt/dr)'^ = c7~^/2 relates the affine parameter to the Lagrangian time. 

C. Geodesic stability 

Using the translation principle, the stability of solutions of natural Lagrangian systems can 
now be discussed in terms of the stability of the corresponding geodesies. It turns out that this 
approach has a number of limitations that are not always appreciated in the literature: 

• the boundaries of the Riemannian manifold; 

• the absence of perturbations that modify energy and the change to affine time; 

• the problem of one-dimensional systems. 

These points equally impact the discussion of global and local stability. 

Boundaries of the Riemannian manifold. The Maupertuis-Jacobi principle translates the so- 
lutions of a natural Lagrangian system (M, L) with given energy E into geodesies of an associated 
n-dimensional metric manifold {M,g^) with metric gE = cr'^k, see (j39p . This Riemannian metric 
has a singular boundary at all points where the conformal factor cr^ vanishes, i.e., for V{x) = E, 
even assuming regularity of the kinetic energy metric k. The Ricci scalar 

= a-^{R^k) - 2{n - 1) □(^j Incr - (n - 2)(n - l)V(fc)'^lna V^^), In^r) , (46) 

generically diverges for a vanishing a. This means that the boundary singularity is a curvature 
singularity, not an artifact due to a possibly poor choice of coordinates. 

The fact that the Riemannian manifold used to model the Lagrangian system has a singular 
boundary can give rise to difficulties in the discussion of the equivalent dynamical system, see 
e.g. 26|. All solutions for which the velocity vanishes at some time parameter will hit this boundary; 



therefore their evolution cannot be determined in the geodesic picture, even though it is perfectly 
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well-defined on the pure dynamical side. Consequently, neither the global nor the local stability 
of solutions hitting the boundary can be determined: the geodesies cannot be extended outside 
a limited range of their proper parameter, so the late time limit that needs to be taken for the 
Lyapunov exponents is ill-defined; also local stability cannot be checked along a complete trajectory. 
Moreover note that fixed point solutions, with vanishing velocity everywhere, are not tractable in 
the geodesic picture at all. This strikingly contrasts the dynamical system analysis, in which 
stability of fixed points is well understood. 




FIG. 2: Potentials of different dynamical systems mapped to the same Riemannian manifold. 

The relevance of these issues is easily illustrated already in simple radially symmetric systems 

L = ]^{x^ + f)-V{r) (47) 

with r = [x^ + y^)^/^. First, consider the harmonic oscillator potential V{r) = and solutions 
of energy E = 1. In this case, the Riemannian manifold will be restricted to {{x,y) | < r < 1} 
with metric g = 2(1 — r'^){dx (iS> dx + dy (i^ dy). The Ricci scalar is given by R = 2(1 — r^)""^, so 
that the boundary is indeed singular. One may check that all radial geodesies, characterized in 
polar coordinates by (/> = 0, reach the boundary in finite proper time. Therefore the manifold is 
also geodesically incomplete. These conclusions hold in full generality for similar potentials, too. 
To illustrate the consequences for the stability analysis we consider the potentials 

V±{r) = 2r2 -r^± 2Q{r'^ - l){r'^ - if , (48) 

where is the Heaviside step function, see figure [2l Solutions with E = \ and initial position 
close to r = are mapped to the same Riemannian manifold for both potentials. So the differences 
of the two dynamical systems cannot be resolved in their geometrical image. In particular, the 
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stability of radial solutions under small perturbations is (under any reasonable intuitive idea of 
stability) dramatically different in the two systems, which the geodesic picture is unable to reveal. 

On the positive side, we remark that in complex systems solutions reaching the boundary of 
the manifold are the exception more than the rule, at least if the kinetic energy metric is positive 
definite. The set of initial conditions for solutions reaching the boundary will be very restricted. 
Solutions with generic initial conditions are expected to be confined to the interior of the manifold. 
In the examples above, only radial geodesies reach the boundary since an arbitrarily small angular 
momentum already prevents points of vanishing velocity along the trajectory. 

Perturbations that modify energy and the change to affine time. The energy of the solutions of 
a Lagrangian system directly enters the definition of the metric of the corresponding Riemannian 
manifold. Hence a generic Lagrangian system cannot be mapped completely into a single metric 
manifold. Note that two metric manifolds {M,gEi) and {M,gE2) corresponding to different choices 
of fixed energy do not possess isometric regions, which makes it impossible to use part of a given 
manifold to discuss dynamics at different energy. For the discussion of both global or local stability 
this means that only perturbations leaving the energy fixed can be considered. From the Lagrangian 
point of view an initial perturbation is represented by a vector ^ in the 2n-dimensional space TpTM, 
where p £ TM provides initial conditions for the reference solution. A perturbation modifies the 
energy of the trajectory unless S,{E) = dE[5) = 0, where E is the energy function defined in 
section llV A[ We can assume dE ^ 0, because the system otherwise would be in a fixed point, so 
dE{^) = restricts ^ to a {2n — l)-dimensional subspace of TpTM. In the geodesic translation we 
therefore expect to have control only over this subspace. 

Another point to note is that the preferred role of the Lagrangian time is lost in translation 
from the Lagrangian picture to the geodesic picture. In the latter the only geometrically natural 
parameter is the affine parameter along the geodesies. That this has impact on the discussion of 
stability becomes clear in simple examples, where time rescalings may render unstable solutions 
stable and vice versa. 

The non-existence of energy-changing perturbations and the loss of Lagrangian time reduce the 
control of stability of natural Lagrangian systems when discussed in the geodesic picture. This fact 
is nicely reflected in the reduction of the number of relevant exponents in the Lyapunov spectrum. 
To see this note that the evolution equation for the perturbation vector in ([9]) can be recast 
as the Jacobi equation (125p using the Levi-Civita connection of gE. In a parallely transported 
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orthonormal basis (cq) = (eo = x, e^) with Ca = the Jacobi equation becomes 

l? = o, ef = -^%/3oef. (49) 

Choosing as initial perturbation = 0, = 0, = gives the solution = const, = 0; 
and the corresponding Lyapunov exponent clearly vanishes. The perturbation corresponding to 
this exponent is an infinitesimal shift along the geodesic. The choice of initial conditions = 0, 
= 0, = gives the solution = const ■ t, (^f = 0, and hence a second vanishing Lyapunov 
exponent. In this case, the initial perturbation corresponds to an infinitesimal shift of the initial 
velocity along the geodesic. As a consequence, only 2n — 2 Lyapunov exponents in the spectrum 
remain significant. 

One- dimensional systems. While it is often argued that the geodesic translation is useful for 

n 

complex systems [8|] , one must emphasize that it completely breaks down for simple one-dimensional 
systems. We already noted above that one obtains two vanishing Lyapunov exponents for any 
translated system; for n = 1 this means that all Lyapunov exponents vanish. So global stability 
cannot distinguish arbitrary one-dimensional systems. The situation is similarly bad for local 
stability. All one-dimensional Riemannian manifolds are flat, so the eigenvalues of the Riemann 
tensor are zero. But then also local stability cannot distinguish arbitrary dynamical systems. Both 
concepts become meaningless for one-dimensional systems. 

After analyzing the different aspects of the geometric picture, we conclude that the Riemannian 
way to define stability suffers from a number of limitations. General dynamical systems cannot be 
mapped into a unique Riemannian manifold, which leads to the presence of singular boundaries 
that may spoil the method of any predictive power. Perturbations that modify the energy cannot 
be discussed at all. Systems with a small number of degrees of freedom are affected more severely; 
the method is not applicable to the one-dimensional case. 

None of these limitations apply to the intrinsic definitions of stability. The generalized Lagrange 
metric is just k without any conformal factor, so that no singular boundaries can occur. Full 
control on perturbations that change energy is expected; indeed, the additional contribution from 
the potential in (j36p in contrast to (09]) removes possible degeneracies that could be responsible 
for a reduced number of Lyapunov exponents. One-dimensional systems do not pose a problem for 
the intrinsic formalism, as we have already seen in the example of the inverted harmonic oscillator. 
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VI. CONCLUSION 

The use of geometric language clearly specifies the intrinsic ingredients of a given physical 
system, and at the same time clarifies which mathematical structures may be meaningfully used. In 
this article we have analyzed the stability of solutions of dynamical systems in geometric language. 
We have restricted ourselves to the important class of continuously evolving systems which can be 
formulated as vector flows X over some smooth manifold S. Although the solutions simply are 
integral curves of the vector field X, this is not a restriction to first order systems; the order rather 
depends on the choice of {S,X). Linear perturbations of reference solutions are characterized by 
vector fields ^ over S that commute with X. Within this basic framework we have reconsidered 
the question of the stability of solutions under linear perturbations. 

For the analysis of global stability exists the extensive mathematical theory of Lyapunov ex- 
ponents. Here the use of the geometric formulation nicely clarifies that global stability in general 
is not an intrinsic notion. It requires outside information that has to be added to the dynamical 
system {S,X). More precisely, a family of seminorms has to be given on the tangent bundle TS. 
This arbitrariness of choice is not sufficiently appreciated in the literature that applies the global 
stability concept. 

In order to analyze the local stability of solutions of (S, X) we have proposed a very intuitive 
new criterion that requires the additional choice of a linear connection on TS to measure local 
accelerations. Hence also local stability in general is not intrinsic to the dynamical system. The 
connection criterion for local stability is applicable to any dynamical system. For the case of second 
order systems we have shown how it reduces to the stability theory of Kosambi, Cartan and Chern. 
In fact, the new criterion is precisely equivalent to KCC theory for natural Lagrangian systems. 

Some dynamical systems contain sufficient geometrical structure to allow for completely intrinsic 
definitions of global and local stability. It is well-known that this is the case for geodesies on a 
Riemannian manifold. But notably, we have shown that this is also the case for general Lagrangian 
systems (M, L). Writing them as dynamical systems (5 = TM,X) shows that the definition of 
global Lyapunov stability proceeds from the generalized Lagrange metric which simply is extended 
to the bundle TS. Local stability according to the newly proposed connection criterion uses the 
Berwald type connection on TS associated to the Lagrangian semispray X. No previous notion of 
local stability seems to cover general Lagrangian systems. 

For natural Lagrangian systems given by Lagrangians quadratic in the velocities, the Jacobi- 
Maupertuis principle is commonly applied to translate solutions of fixed energy into geodesies of 
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an associated Riemannian manifold. We have compared the resulting notions of geodesic stability 
with the intrinsic notions derived from the Lagrangian system. We have pointed out that geodesic 
stability suffers from various illnesses: boundaries of the Riemannian manifold that are curvature 
singular; the non-ability to resolve perturbations that modify energy; and the complete breakdown 
of the analysis for one-dimensional systems. None of these issues troubles the intrinsic definitions 
of stability in the Lagrangian picture. 

The main achievement of this article lies in clarifying in a geometric framework the structures 
that underlie any discussion of the global or local stability of solutions of dynamical systems. It 
has emerged that the geometry of dynamical systems forces us to define a new criterion for local 
stability in terms of a choice of linear connection. Our comparison of this criterion to the geodesic 
stability of natural Lagrangian systems already shows a number of advantages of the intrinsic 
geometric approach. We expect that the practical relevance of our construction will be proven in 
future applications. 
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